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We study analytically the spacelike singularity inside a spherically-symmetric, charged black hole 
coupled to a self-gravitating spherical massless scalar field. We assume spatial homogeneity, and 
find a generic solution in terms of a formal series expansion. This solution is tested against fully- 
nonlinear and inhomogeneous numerical simulations. We find full compliance between our analyt- 
ical solution and the pointwise behavior of the singularity in the numerical simulations. This is a 
strong scalar-curvature monotonic spacelike singularity, which connects to a weak null singularity 
at asymptotically-late advanced time. 



PACS number(s): 04.70.Bw, 04.20.Dw 
I. INTRODUCTION 

The spacetime singularities, which classical General 
Relativity (GR) predicts under very plausible assump- 
tions to inevitably occur inside black holes |Q] , are a ma- 
jor challenge for Physics, as the currently known physical 
laws are presumably invalid at such singularities. In- 
stead, some other theory, as yet unknown, is expected to 
take over and control their structure. However, the GR 
predictions are still of the greatest importance, as they 
reveal the spacetime structure under extreme conditions 
in the strong-field regime. 

In the last few years, much advance in the study of 
black hole singularities has been achieved in the frame- 
work of classical GR. Poisson and Israel showed that 
when non-linear perturbations were allowed, the inner 
horizon of Reissner-Nordstrom evolved into a scalar cur- 
vature singularity, where the internal mass function di- 
verged and curvature blew up For spinning or 
charged black holes, it has been found that at least for 
a portion of the singularity the latter is null and weak 
HI). Namely, despite the divergence of spacetime cur- 
vature at the singularity, its tidal influence on extended 
physical objects is bounded. The detailed structure of the 
null weak singularity was studied numerically and 
analytically in the context of a spherical charged 
black hole, non-linearly perturbed by a self-gravitating, 
minimally-coupled, massless scalar field. (In the case of 
a spinning black hole, the evidence for the occurrence of 
the null weak singularity emerged from analytical per- 
turbative [Q and non-perturbative || analyses.) It has 
been found numerically H§], that under the impact of 
the non-linear scalar field, the generators of the Cauchy 
horizon are monotonically focused, until the focosing ef- 
fect is complete, and the singularity becomes spacelike, 
rather than null. 

Although the existence of a spacelike singularity was 
suggested by several independent numerical simulations 
jlO|,|]]|] , that spacelike singularity has not been shown 
analytically to exist: its very occurrence is a non-linear 



effect, and no fully non- linear analytical analyses of the 
problem have been performed. In addition, the proper- 
ties of that singularity (both geometrical and physical) 
have not been studied in detail. In fact, the occurrence 
suggested by the numerical simulations is the only thing 
currently known about the spacelike singularity inside 
charged black holes. It is the purpose of this paper to 
study the spacelike singularity in the model of a spheri- 
cal charged black hole and a neutral massless scalar field 
(the same model for which the singularity was found nu- 
merically). 

The organization of the paper is as follows. In Section 
H we show that within a simplified homogeneous model 
one indeed finds a generic solution describing a spacelike 
singularity. This model is used to analyze the behavior of 
geometry and the scalar field near the singularity. In Sec- 
tion III we describe fully-nonlinear and inhomogeneous 



numerical simulations which test the validity of the ho- 
mogeneity assumption of our analytical model. We find 
full compliance between the analytical results and the 
pointwise behavior found in the numerical simulations. 
We also use the numerical simulations to study the spa- 
tial structure of the singularity. We conclude and discuss 



our results in Section IV 



II. ANALYTICAL HOMOGENEOUS MODEL 

In order to facilitate the analysis, we make the sim- 
plifying assumption, that the singularity can be mod- 
eled to be homogeneous. (This assumption is justified a 
posteriori numerically.) That is, we assume that spatial 
gradients of the dynamical fields are negligible near the 
spacelike singularity compared with temporal gradients. 
Consequently, our analysis can actually probe just the 
pointwise behavior at the singularity. The spatial depen- 
dence, however, cannot be studied with this model. We 
shall therefore study the spatial dependence numerically 
below in Sec. III. 
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We write the general homogeneous spherically- 
symmetric line element as 

ds 2 = h{r)dt 2 + f{r)dr 2 + r 2 dtt 2 (1) 

where dVl 2 = d6 2 + sin 2 9 dip 2 is the metric on the unit 
two-sphere. Here, r is the radial coordinate, defined such 
that spheres of radius r have surface area Airr 2 (note 
that r is timelike inside the black hole), and t is nor- 
mal to r (t is spacelike inside the black hole. There is a 
gauge freedom in t — see below). The Einstein-Maxwell- 
scalar equations (with a free electric field corresponding 
to charge q and with a scalar field <&) are then given by 
(/V + / 2 -/)/(/r 2 )=cf>'2 + /(? 2 /r 4 (2) 



(h'r - hf + h)/ (hr 2 ) = $' 2 - fq 2 /r 4 
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I (Arfh 2 ) = fq 2 y - $' 2 (4) 

in addition to the Klein-Gordon equation VV^ $ = 0, 
whose first integral is $'(r) = d y/\f] (hr 4 ) |, where d ^ 
is an integration constant. Here, a prime denotes differ- 
entiation with respect to r, and V M denotes covariant 
differentiation. We seek a generic solution to these equa- 
tions describing a spacelike singularity at r = 0. 

The source terms for the field equations [the right-hand 
side (RHS) of Eqs. (0)-(Q)] contain contributions from 
both the electric field and the scalar field. We find that 
for a dominant electric-field contribution near the singu- 
larity at r — 0, the latter is timelike, rather than space- 
like. As in this case the scalar field's contribution to the 
field equations is negligible compared with the electric 
field's contribution near r = 0, to the leading order in 
r one would expect to find the same solution as with no 
scalar field at all. However, from the generalized Birkhoff 
theorem the solution in the latter case is nothing but 
the Reissner-Nordstrom solution. Namely, with domi- 
nant electric field, to the leading order in r the singu- 
larity is indistinguishable from the Reissner-Nordstrom 
singularity, which is timelike. We note in passing that 
this case might be of some relevance for consideration of 
the self-gravitating scalar field for a hypothetical exten- 
sion of the spacetime manifold beyond the weakly singular 
Cauchy horizon. A second case is scalar field and elec- 
tric field which are comparable in strength near r = 0. It 
turns out that there is no consistent solution for compara- 
ble electric- and scalar-field contributions at the singular- 
ity^. The third case is the case where the electric field's 



*We remark that this does not mean that there cannot be a 
solution with comparable contributions of both fields. We do 
find, however, that there is no homogeneous singularity with 
comparable contributions. 



contribution to the energy-momentum tensor near the 
singularity is negligible compared with the scalar field's 
contribution, and to this case we shall specialize below. 
We emphasize that the causal structure of the singularity 
at r = depends on which field dominates: if the elec- 
tric field dominates the singularity is timelike, and if the 
scalar field dominates the singularity is spacelike. Note, 
that in this sense the spacelike singularity we are study- 
ing is similar to the spacelike singularity with vanishing 
charge. In fact, we find that to the leading order in r the 
singularity in our case is indistinguishable from the sin- 
gularity in the uncharged case. However, in our case the 
spacelike singularity is a non- linear effect, whereas in the 
uncharged case a spacelike singularity exists even with 
a vanishing scalar field (the Schwarzschild singularity). 
(There are additional differences between the spacelike 
singularities in the two cases pd| .) Another comment 
to be made about the timelike singularity is that it is 
reasonable to expect a linear perturbation analysis to be 
valid near it, as the non- linear effects of the scalar field 
are negligible. 

We assume a formal series expansion for the metric 
functions of the general form f^ n \r) — Y^ii=i /«( r ) ano - 
h^ n \r) = ^i( r )j an d assume that both series have 

a finite radius of absolute convergence to / and h, re- 
spectively, in the limit n — > oo. We also assume that for 
all values of i, fi(r) = fior mi and hi(r) — h i0 r ni , with 
constant fio and /i^o- (These assumptions are justified a 
posteriori numerically.) We then find the leading terms 
in r for the metric functions, namely 

+ l) 2 (3/5 + 2) 
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h(r) = d 2 C r? - d 2 {P + l M + 2) C 2 q 2 r 2 ' 3 + 0(r 2 ^) 

P 

(6) 

where C and [3 are positive parameters. Here, 7 = 2 if 
(3 > 2 and 7 = [3 otherwise. For the scalar field we find 

(/3+l) 3 / 2 (3/3 + 2) c2 q2 r0 



$(r) = V/3+1 hir 
+ 0{r^). 



2/3 3 



(7) 



This form for the solution suggests that /, h are given by 
infinite double series 
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f,h=j2E a ™ rf3m+(0+2){n ~ 1)+ai ' h > ( 8 ) 

m— 1 n— 1 

aLn = a mn(C>/3) and a h mn = d 2 a ] } nri (C, (3) being the ex- 
pansion coefficients. Here, <jj — and ah = —2. [|] As 



'We note that whereas with vanishing electric charge it has 
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we have already remarked, the black hole electric charge 
does not appear in the leading order terms of either the 
metric functions or the scalar field. This is a direct conse- 
quence of the dominance of the scalar field over the elec- 
tric field near the singularity. Consequently, the deriva- 
tion of the solution is similar to the derivation in , and 
further details can be found there. The charge affects the 
spacetime geometry and the behavior of the scalar field 
in higher order terms, however, as is evident from the so- 
lution, such that at any finite value of r, even very close 
to the singularity, both the geometry and the field are 
different from the uncharged case's counterparts. Note, 
that we find the parameter (3 to be positive. In the case of 
the spacelike singularity inside an uncharged black hole, 
the corresponding parameter is found to be greater than 
— 1, which is the value found for the Schwarzschild sin- 
gularity. Thus, in the uncharged case the parameter (3 
changes continuously from its vacuum value, whereas in 
the charged case we study here it jumps discontinuously 
from the —2 value of the Reissner-Nordstrom singularity 
which corresponds to the (electro-)vacuum limit of our 
spacetime, to positive values. This discontinuity is the 
result of the non-linear nature of the spacelike singular- 
ity: with no scalar field (the electro-vacuum case) the 
singularity at r — is timelike. With the addition of a 
scalar field the only singularity we find is spacelike. The 
change from timelike to spacelike is discontinuous, and 
this is what we indeed find. 

From this solution, we draw the following inferences: 
First, we define a 'tortoise' coordinate by g r , r , = —gtt- 
Then, null coordinates are defined by t — (v — u) /2 and 
r* = (v + u)/2, with ingoing v and outgoing u. For a 
certain outgoing ray u = uq = const which runs into 
the spacelike singularity, the latter is hit at some finite 
value of advanced time v = (see Fig. |l|). (The null 
singularity is located at v — oo.) We find that to the 
leading order in — v 



spacelike singularity 



r(u ,v)= [d 2 /(p + l)] 1/4 („.-«) 



a/2 



(9) 



This property of the spacelike singularity inside a non- 
linearly perturbed charged black hole is similar to the 
properties of the singularity inside uncharged black holes 
|pl| . We next consider an observer who follows a radial 
t = const trajectory. Let r be that observer's proper 
time, set such that r = at r = 0. We denote by x° the 



been shown that the expansion coefficients in the series ex- 
pansion analogous to (^) can be found uniquely for all orders 
[ p^[ , in the case studied here this is yet to be proved. However, 
we believe that this is also the case here, despite the formal 
difference in the series expansions of the two cases: whereas 
in the uncharged case the expansion is given in the form of 
a simple series, in the charged case it is given in terms of a 
double series. 
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FIG. 1. The Penrose diagram of the simulated spacetime. 
Singularities are depicted by thick lines. Dashed thick lines 
denote regions of the singularity not covered by the domain of 
integration, and describe a possible picture of these portions. 
The causal structure of the singularity changes from null to 
spacelike at (u = u, v = oo). For any outgoing null ray at 
u = no > u the spacelike singularity is reached at v = «*. 



coordinate tangent to the worldline. Then, to the leading 
order in r, the (00) tetrad component of the Ricci tensor 
and the Kretschmann scalar, correspondingly, are given 
by 



r (o)(o)(t) = 8 



R aPl sR a ^ S {T) = 64 



(3+1 1 
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(10) 



(11) 



As Eqs. ([!(]) and (^lj) are given to the leading order 
in r, they are identical to the corresponding equations 
in the uncharged case (cf. |ll]]). Equation (11) implies 
that this is a scalar-curvature singularity. From Eq. ( p^| ) 
and from the theorem by Clarke and Krolak [jl3| it then 
follows that this singularity is strong in the Tipler sense 
ijlij , as the twice integrated R(o)(o)( T ) over T diverges 
logarithmically as r — > 0. Namely, any extended physical 
object will unavoidably be crushed to zero volume upon 
arrival to the singularity. As we find that (3 > 0, we 
infer that there is no finite (3 which nullifies any of the 
RHS's of either Eq. © or Eq. (0). Consequently, 
for the entire range of permissible parameters, this is a 
strong scalar-curvature singularity. This property of the 
singularity is in sharp contrast with the Cauchy horizon 
singularity, which is known to be null and weak. 

Let us consider the genericity of our solution. The 
notion of a general solution for non-linear differential 
equations is not unambiguous. However, we can count 
the number of arbitrary parameters in our solution, and 
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thus show that it is generic. From the physical point of 
view, one would expect two arbitrary parameters (one 
because of Birkhoff 's theorem and one due to the scalar 
field. Note that as the scalar field is neutral, the electric 
charge is fixed). Apparently, our solution has three ar- 
bitrary parameter, namely /?, d 2 , and C. However, the 
arbitrariness in d 2 reflects merely a trivial gauge mode, 
i.e., the freedom to re-scale the coordinate t (see Ref. 
jnj). Consequently, our solution has two physical de- 
grees of freedom, and is therefore generic. These consid- 
erations do not rule out the possible existence of other 
generic solutions. However, this solution is the one which 
is realized in the numerical simulations (see below). 

III. FULLY INHOMOGENEOUS NUMERICAL 
SIMULATIONS 

The analytical approach presented above assumes the 
singularity to be homogeneous. However, the spacetime 
inside black holes (even very close to the singularity) is 
expected in general to be inhomogeneous. To what extent 
then is the homogeneity assumption restrictive? We shall 
present in what follows the results of fully non-linear and 
inhomogeneous numerical simulations. We show, that 
the above homogeneous analytical model succeeds to de- 
scribe the pointwise behavior of the singularity as inferred 
from the inhomogeneous simulations. In this sense the as- 
sumption of homogeneity is a simplification which does 
not ruin the structure of the singularity, as the numeri- 
cal simulations reveal that the dependence of the geom- 
etry at and near the singularity depends only weakly on 
spatial coordinates: it is the temporal gradients which 
dominate. However, the spatial dependence can still be 
studied numerically. We shall do just this after we con- 
front the predictions of the analytical model with the 
pointwise behavior of the numerical simulations. 

Let us discuss then the results obtained from fully non- 
linear (and inhomogeneous) numerical simulations. We 
used the same code which was used in Ref. [|| (see also 
flBf ) to study the null singularity in this model. This code 
is based on double-null coordinates and on free evolution 
of the metric functions and fields. The code is stable 
and converges with second order. (Slight adaptations of 
the code were needed, to allow for a careful approach to 
r = 0. These modifications are described in Ref. [prf .) 
Our initial value set-up is such that prior to the smooth 
initial pulse of compact support the geometry is Reissner- 
Nordstrom with charge q and mass M (see Fig. Q). Then, 
due to the scalar field the geometry is changed, and the 
final mass of the black hole is Mf. 

Figures ||, || and [| display the results of the numer- 
ical simulations with q/M = 0.998 and M f /M w 1.51. 
We have also performed numerical simulations with other 
values of the parameters, and found no qualitative change 
of the behavior. In fact, when we increase the amplitude 
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FIG. 2. (A) The metric function G = -2g uv . (B) The 
scalar field $. (C) The parameter (3 as calculated from the 
exponent of h. (D) The ratio of (3 as calculated from the 
exponent of h and j3 as calculated from the amplitude of 
In all graphs the abscissae are r. 

of the scalar-field pulse (while keeping the support of the 
pulse fixed), we find the following two effects. First, the 
ratio Mf/M increases, as should be expected from the 
greater energy content of the pulse. Second, the larger 
the amplitude, the earlier the spacelike singularity oc- 
curs. Namely, with larger amplitude for the initial pulse 
the generators of the Cauchy horizon focus faster, and 
the Cauchy horizon focuses completely within a smaller 
lapse of affine parameter. However, the spacelike sin- 
gularity itself is unchanged by this change of the initial 
parameters, and is in this sense robust. Very similarly, 
varying the value of q/M does not change the properties 
of the singularity. This should indeed be expected from 
our analytical considerations of Section [n| the singular- 
ity is dominated by the scalar field, and to the leading 
order in r the charge of the black hole has no effect. 

The data are shown along an outgoing ray at Uq ap- 
proaching the spacelike singularity at r = (similar re- 
sults are obtained for any choice of such an outgoing ray; 
in all four graphs of Fig. || the abscissae are r): §(A) 
the metric component g uv , for double- null coordinates 
defined such that v = r on the outgoing segment of the 
characteristic hypersurface, and u cx — r on its ingoing 
segment. The vanishing of g uv implies the vanishing of 
both / and h. Moreover, the asymptotically linear curve 
in this plot implies a power-law behavior near the singu- 
larity. (This is checked more accurately in Fig. |](C).) 
Figure |](B) displays the scalar field <F We find that $ di- 
verges logarithmically, in agreement with Eq. (Q) . Figure 
§(C) shows the local power index of g uv , namely, the lo- 
cal slope in Fig. §(A). The asymptotically constant value 
to which the local power approaches — in agreement with 
Eqs. (||) and (||) — implies that the power-law assumption 
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FIG. 3. The radial coordinate r(uo,v) as function of (ad- 
vanced-time) distance v* — v from the spacelike singularity. 
The singularity is located at v = v, . 



200 
180 
160 
140 
120 
S 100 

-Q 

80 
60 
40 



20 



-_ 

6.5 



CP" 



FIG. 4. The value of /3 for different outgoing null rays as a 
function of w* , namely, as a function of the value of the ingoing 
coordinate v for which the spacelike singularity is reached 
when approached along outgoing null rays. 



for the metric functions is justified. In fact, this power- 
law index is just the value of the parameter (3 for the in- 
tersection point of the outgoing ray under consideration 
with the spacelike singularity. Finally, we test another 
aspect of the analytical model. We found that the same 
parameter (3 appears in two different contexts, namely, in 
the exponent of the metric functions [Eqs. (|J) and (ph], 
and in the amplitude of the scalar field [Eq. (0)]. The 
ratio of these two values for (3 is predicted, therefore, to 
asymptotically approach a value of unity. This is indeed 
shown in Fig. ||(D), to within 1 part in 10 3 . 

Figure [| shows r(uo, v) as a function of the distance (in 
terms of advanced time) from the singularity — v. The 
slope is equal asymptotically to the value of 1/2 implied 
by Eq. @ to within a numerical error of 1%. 

Finally, we use the numerical simulation in order to 
check the spatial dependence of the solution. Figure |] 
displays the value of the parameter j3 as a function of 
advanced time, namely, as a function of the value of u* 
at which the outgoing ray under consideration hits the 
singularity. We find that along the spacelikc singularity, 
upon approach to the spacetime event where the causal 
structure of the singularity changes from null to spacelike, 
(3 increases rapidly. This suggests that at the particular 
outgoing ray at u = u for which that event is reached [3 
may be infinite. We can therefore compare the relative 
strength of different points along the spacelike singular- 
ity. Although this singularity is strong throughout, the 
nearer one is to the spacetime event where the causal 
structure changes, the weaker is the strong singularity, as 
curvature scalars diverge slowlier. That is, let us consider 
a constant r hypersurface near the singularity (namely, 
deep in the asymptotic region). Then, approaching u 
along this hypersurface curvature scalars decrease (as (3 
increases). This fits well with our overall picture: when 



the singularity first forms it is null and weak. At later 
times, the weak null singularity becomes stronger with 
increasing affine parameter. This is manifested by the 
focusing effect, and by the inverse proportion of the blue- 
shift factors to the area coordinate ||. Thus, we have 
a weak singularity which becomes stronger. Then, at 
u — u, r vanishes and the causal structure of the singu- 
larity changes and it becomes spacelike and strong, and 
strengthens as one gets farther from the changing point, 
namely, as u — u increases. 

IV. DISCUSSION 

We find analytically a simple homogeneous model, 
which describes the pointwise behavior of the spacelike 
singularity which occurs inside spherical charged black 
holes with a non-linear scalar field. This singularity is 
generic in the sense that the solution has the expected 
number of free parameters. The numerical simulations 
imply that the singularity we find in the simple homoge- 
neous model successfully describes the pointwise behav- 
ior of the singularity of the fully nonlinear and in general 
inhomogeneous spacetime. The spacelike singularity is 
scalar curvature (like the null singularity which precedes 
it), monotonic, and strong (unlike the null singularity, 
which is weak). 

It remains an open question, to what extent this simple 
model captures the essence of the singularity inside realis- 
tic black holes, namely, spinning black holes with vacuum 
perturbations: First, scalar fields are known to be related 
to unique phenomena, such as the destruction of the oscil- 
lations of the Belinskii-Khalatnikov-Lifshitz (BKL) sin- 
gularity p6| . Second, we find that the scalar field in our 



5 



model is dominant near the singularity, whereas the BKL 
singularity is often dominated by vacuum perturbations 
rather than by matter fields . Finally, both the spher- 
ical symmetry and the scalar field are mere toy models 
for the more realistic cases. A major open question is 
then whether in more realistic cases there would indeed 
be a spacelike singularity to the future of the null sin- 
gularity, and whether it would be monotonic - like the 
spacelike singularity inside spherical charged black holes 
with a scalar field - or oscillatory - like the BKL singu- 
larity. 
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